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Abstract 

We carried out two-dimensional magnetohydrodynamic simulations of the Galactic 
gas disk to show that the dense loop-like structures discovered by the Galactic center 
molecular cloud survey by NANTEN 4 m telescope can be formed by the buoyant 
rise of magnetic loops due to the Parker instability. At the initial state, we assumed 
a gravitationally stratified disk consisting of the cool layer (T ~ 10 3 K), warm layer 
(T ~ 10 4 K), and hot layer (T ~ 10 5 K). Simulation box is a local part of the disk 
containing the equatorial plane. The gravitational field is approximated by that of 
a point mass at the galactic center. The self-gravity, and the effects of the galactic 
rotation are ignored. Numerical results indicate that the length of the magnetic loops 
emerging from the disk is determined by the scale height of the hot layer (~ 100 pc 
at 1 kpc from the Galactic center). The loop length, velocity gradient along the loops 
and large velocity dispersions at their foot points are consistent with the NANTEN 
observations. We also show that the loops become top-heavy when the curvature of 
the loop is sufficiently small, so that the rising loop accumulates the overlying gas 



faster than sliding it down along the loop. This mechanism is similar to that of the 
formation of solar chromospheric arch filaments. The molecular loops emerge from 
the low temperature layer just like the dark filaments observed in the Ha image of 
the emerging flux region of the sun. 

Key words: Galaxy: Magnetic loops, ISM - - Sun: Chromosphere, Magnetic 
fields, Magnetohydrodynamics 

1. Introduction 

Spiral galaxies have large-scale mean magnetic fields. It is widely accepted that the 
magnetic fields are amplified and maintained by the dynamo mechanism (e.g, Parker 1971). 
Mathewson and Ford (1970) showed by observations of optical polarization due to dust grains 
that the magnetic fields are nearly parallel to the galactic plane with wave-like patterns or loop- 
like structures, suggesting the emergence of magnetic field via magnetic buoyancy (e.g., Parker 
1966; Mouschovias et al. 1974; Blitz and Shu 1980; see also Tajima & Shibata 1997). A typical 
magnetic field strength in the Galactic plane is a few \i G. The magnetic field strength increases 
with the increase in density toward the Galactic center. According to the radio observations 
of the Galactic center and its vicinity, the magnetic field strength reaches a few mG in local 
regions near the Galactic center (Morris & Serabyn 1996). 

The Galactic gas disk consists of a gas with three phases, namely, cold (~ 30 K), warm 
(~8x 10 3 K), and hot (~ 10 6 K) components (Mckee & Ostriker 1977). The cold component 
corresponds to the molecular clouds observed with radio molecular lines (such as CO lines), 
while the warm component is observed with the Hi line, and the hot component is observed 
with soft X-rays. 

The Galactic molecular cloud survey by NANTEN 4 m telescope (Mizuno & Fukui 2004) 
found two dense gas features having a loop-like shape with a length of several hundred pc and 
width of ~30 pc within 1 kpc from the Galactic center (Fukui et al. 2006). Fukui et al. (2006) 
discovered large gradients in the line-of-sight velocity (~ 30 km s" 1 ) along the molecular loops. 
Moreover, large velocity dispersion (~ 30 — 50 km s _1 ) is observed near the foot-points of the 
loops. Since the loops have a total mass of ~ 1O 5 M , the kinetic energy of a loop is estimated 
to be ~ 10 51 erg. 

This energy is too large to be explained by a single supernova explosion. Moreover the 
velocity distribution is distinct from that of an expanding shell. These features can be explained 
by the magnetic loops buoyantly rising due to the Parker instability (Parker 1966). 

Matsumoto et al. (1988) carried out two-dimensional (2D) magnetohydrodynamic 
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(MHD) simulations of the Parker instability. They found that dense regions are formed in 
the valley of magnetic loops in the nonlinear stage of the Parker instability and that shock 
waves are formed at the footpoints of the rising magnetic loops where the gas infalling along 
the magnetic field lines collides with the dense gas near the equatorial plane of the disk. Shibata 
& Matsumoto (1991) applied the results of the MHD simulations to the formation of molecular 
clouds in the Orion region. Subsequently, Kamaya et al. (1996) studied the triggering of the 
Parker instability by supernova explosions. Basu, Mouschovias & Paleologou (1997) investi- 
gated the effect of the Parker instability on the structure of the interstellar medium, and Kim 
et al. (2000) applied it to galactic disks. Machida, Hayashi & Matsumoto (2000) reported 
the results of the three-dimensional (3D) global MHD simulations of the Parker instability in 
a differentially-rotating disk. They found that magnetic loops emerging from the disk form a 
structured corona above the disk. 

Fukui et al. (2006) presented the results of 2D MHD simulations of local volume of the 
Galactic disk. They suggested that the downflows can be the origin of the violent motion 1 and 
extensive heating of the molecular gas in the Galactic center. In their simulations, however, 
the cold component of the galactic gas that corresponds to molecular clouds was not taken into 
account and the mechanism of the formation of dense loop-like structures was not clear. 

In the solar atmosphere, arch filaments are observed in the Ha images of the chromo- 
sphere. Isobe et al. (2006) found that dense filaments similar to Ha arch filaments are formed 
in the emerging flux. These filamentary structures are cool and dense above the chromosphere. 
This configuration in which the dense regions exist around the top of the emerging loops is a 
typical feature of the emerging loops and is a good approximation to what will be observed 
as dark features in the Ha image. They clarified the reason why the emerging loop becomes 
top-heavy on the basis of the results of 3D and 2D MHD simulations of the emergence of the 
magnetic loops from the convection zone, through the chromosphere, to the corona. The dense 
gas accumulated around the top of the rising loops fragments into filaments by the Rayleigh- 
Taylor instability, and slides down along the magnetic field lines. These filaments correspond 
to the arch filaments observed in the emerging flux regions of the sun. 

In this paper, we present the results of 2D MHD simulations of the Galactic gas disk. 
We will consider a local volume of the Galactic gas disk that consists of a low temperature 
layer (the cold component of the Galactic gas disk), warm layer and hot layer through 2D 
MHD simulations. In section 2, we will present basic equations and the numerical model used 
in this study. The numerical results will be reported in section 3 and section 4 is devoted for 
discussion. Summary will be given in section 5. 

1 Among the mysterious features unique to the Galactic center, the violent motion in the central molecular 
zone and the other similar broad velocity features whose velocity dispersion is about 15 to 30 km s _1 have 
been highly enigmatic for the past few decades. 
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2. Numerical Model 



2.1. Assumptions and Basic Equations 

We adopt the local Cartesian coordinates (x,y,z) at radius Tq from the Galactic center 
(see Figure 1). The x-direction is taken to be the azimuthal direction, and the z-direction is 
parallel to the rotational axis of the Galactic disk. The partial derivatives of the background 
medium with respect to the radial direction y are neglected (i.e., local approximation). The 
assumptions, basic equations and initial conditions are similar to those in Matsumoto et al. 
(1988). We will assume the following: (1) the medium is an ideal gas with the specific heat 
ratio 7 = 1.05, (2) the viscosity and resistivity are neglected, (3) the effects of the rotation of 
the disk and self-gravity are neglected. 

The basic equations are: 

^ + (v- V)p=-pV-v, (1) 

dv 11 

+ (v-V)v = —VP + — (V x B) x B + g, (2) 

at p Anp 



and 



OP 
~dt 

dB 



(«.V)P=-7P(V-t>), (3) 

m Vx(nB). (4) 

Here, g = (0,0, g(z)) is the gravitational acceleration assumed to be 

g{z) = -GMz/{rl + z 2 fl 2 , (5) 

where G is the gravitational constant. This vertical component of the gravitational acceleration 
is produced by the gravitational potential by a point mass and gives a simple but exact model 
for the accretion disk rotating around the point mass. In the case of galaxies, this is not an 
exact model, but at least the region of the equator can be approximated by this model, because 
the gravitational acceleration near the disk in galaxies is approximately proportional to the 
height from the equatorial plane of the disk. The radial component of gravity is assumed to be 
equal to the centrifugal force due to the rotation of the disk. The effects of the rotation will be 
discussed in section 4.4. 

Equations (l)-(4) are rendered dimension less by using normalizing constants r , Cso 
and po, where Cso is the sound speed in the mid-temperature region and po is the unperturbed 
density at the equatorial plane. When the numerical results are compared with observations, 
we use the units of length, velocity, and time in the simulation to be r = 1 kpc, Cso — 18 km s _1 , 
and to = r / Cso = 5.6 x 10 7 yr, respectively. The unit temperature is T = pC 2 Q / (jR) = 2 x 10 4 K, 
where p and R are the mean molecular weight and gas constant, respectively. The normalization 
units are summarized in table 1. 

We also introduce a non-dimensional parameter e = 7(V^/Cf ), where Vk = (GM/ro) 1 ^ 2 
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is the Keplerian velocity at radius tq. In this simulation, we adopt e = 130 and Vk = 200 km s _1 . 
For this parameter, the sound speed in the low-temperature region is about 5.6 km s" 1 whose 
value is close to the observed velocity dispersion (5 — 9 km s" 1 ) of the main gas components 
(Boulares & Cox 1990). However, this sound speed is larger than that of the molecular gas. 

2.2. Initial State 

The initial state is assumed to be in magnetohydrostatic equilibrium. The gas layer is 
initially composed of three layers: the cool equatorial layer (T = T c , \z\ < Zi), the warm (mid- 



temperature) layer (T = T m ,z\ < \z\ < z 2 ), and the hot galactic halo (T 
initial distribution of temperature is assumed to be 

\z\ — Z\ 



T h , \z\ > z 2 ). The 



T(z)=T c + (T n 



+ (T h 




tanh 




'"'^ ~ z ]+!}•■ (6) 

In this study, we take T c = 0.1T , T m = T , Th = 10T , w t = 0.015r , z\ = 0.12r and z 2 = 0.24ro- 
Although these values are not realistic for the Galactic gas disk, they would be acceptable 
for our first attempt to study the effects of the cold component of the galactic gas and the 
mechanism of the formation of dense loop-like structures. Numerical results do not depend 
much on the temperature of the hot galactic halo (Kamaya et al. 1997). Further discussion on 
the dependence of numerical results on the disk temperature will be given in section 4.5. 

We assume that the magnetic field is initially parallel to the equatorial plane; B = 
(Z? x (z),0,0), and is localized in the cool equatorial region (\z\ < z\). The magnetic field strength 
is determined by introducing the plasma beta (the ratio of gas pressure to magnetic pressure, 
hereafter denoted by (3) as 



BJz] 



]txP(z) 



\ 00 



where 




tanh 



Zi 



+ 1 



(7) 



(8) 



Here, (3q is the plasma (3 at the galactic plane, and Zf is the half thickness of the magnetic flux 
sheet. The initial density and gas pressure distributions are calculated numerically by solving 
the equation of magnetohydrostatic equilibrium: 



d 

dz 



P(z) 



Bj{z) 
8tt 



+ p{z)g{z) = Q. 



(9) 



We have studied the evolution for five models, whose parameters are summarized in Table 
2. The distributions of the initial temperature T, density p, gas pressure P, and magnetic 
pressure B 2 x j (8ir) are shown in Figure 2a. Figure 2b shows the profile of the initial local pressure 
scale height A(z) = Cs{z) 2 /( / -fg(z)) and the pressure scale height including the magnetic field 
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Ab(-z) = (1 + /5(^) _1 )A(z) for model B ((3 = 1), which we call the fiducial model hereafter. 
2. 3. Stability of the Equilibrium Model 

The equilibrium state we described in section 2.2 is unstable against the Parker insta- 
bility. The Parker instability in an isothermal gas has the linear maximum growth rate at a 
finite wavelength Ap ar ker ~ 10A — 20A, where A is the scale height (Parker 1966). This is be- 
cause small- wavelength modes are stabilized by the magnetic tension force. We analyzed the 
linear stability of the initial model with a normal mode method similar to that of Horiuchi 
et al. (1988). The linearized equations are the same as those in Horiuchi et al. (1988) and 
Kamaya et al. (1997). We consider the growth of a small perturbation that has a functional 
form 5W oc exp(iut + ik x x + ik y y), where W is a physical quantity (p,P,v z ,B x ) and SW is its 
perturbation. The eigenvalues (u) and the eigenfunctions are calculated numerically. 

Figure 3a shows the linear growth rate (iu) of the fundamental mode (Horiuchi et al. 
1988) of the Parker instability as a function of the horizontal wave number k x for five cases, 
/3q — 0-5, 1,2,4 and 10 when Zf/r = 0.08 and k y = 0. Figure 3b shows that the maximum growth 

— 1/2 

rate is inversely proportional to the square root of (3 (iu oc /3 ). This result is consistent 
with that of the linear analysis of the Parker instability in uniform gravitational fields (Parker 
1966). The most unstable wavelength (A max ) is A max — (7r/8)r for model B ({3 = 1)- It is noted 
from Figure 2b and Figure 3 that A max is nearly 10 times the local pressure scale height of the 
mid-temperature region in 0.12 < z/r < 0.24 (A max ~ 10A ~ 400 pc) and is much larger than 
that of the low-temperature cool region around z/ro ~ 0.1. 

2.4- Boundary Conditions and Numerical Method 

We assumed free boundaries at z — Z m \ n and z = Z max such that waves transmit freely 
by setting ^-derivatives of all the variables vanish, and imposed periodic boundary conditions 
at x = ATmin and x = X max . 

In order to trigger the Parker instability, small velocity perturbations of the form 



are given initially within a finite horizontal domain (\x\ < A x /2) on the magnetic flux sheet, 
where the perturbation wavelength X x is close to that of the most unstable wavelength for the 
Parker instability. Here A (= 0.01) is the maximum value of v z /Cso in the initial perturbation. 

The numerical scheme we used is the Rational CIP (Cubic interpolated profile) method 
(Yabe & Aoki 1991; Xiao, Yabe & Ito 1996) combined with the MOC-CT method (Evans & 
Hawley 1988; Stone & Norman 1992). The magnetic induction equation was solved by the 
MOC-CT and the other equations were solved by the CIP (Kudoh, Matsumoto & Shibata 
1998, 1999). 

The size of the simulation box is (A" max — X m[n , Z max — Z min ) = (4r , 4r ) = (19 H, 19H), 
where H is the scale height at the point where the gravity is maximum (H = Cg/(7g max ) = 




(10) 
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0.21r ). The grid sizes are Aa; = 2.5 x 10 3 r for |a;| < 0.8r , Az = 2.5 x 10 3 r for \z\ < 0.8r 
and they slowly increase for |x| > 0.8ro, or \z\ > 0.8ro by an increment of 5% at each grid (e.g, 
|Ax m | = 1.05|Ax;|). The number of grid points is (N X ,N Z ) = (472,472). 

3. Numerical Results 

Figure 4 shows the time evolution of the fiducial model. Solid curves depict magnetic 
field lines. Color shows the density distribution. Arrows show velocity vectors. The overall 
evolution agrees with that of Matsumoto et al. (1988). That is to say, as the instability grows, 
the magnetic field lines bend across the equatorial plane. As the gas slides down along the 
undulating magnetic field lines, the rarefied regions buoyantly rise, and form magnetic loops in 
the later (nonlinear) phase. In the valleys of the magnetic loops, dense spur-like structures are 
created almost perpendicular to the galactic plane. At the top of the emerging loops, dense 
shell-like structures are formed (Figure 5a). These structures were not recognized in previous 
simulations (e.g., Matsumoto et al. 1988) partly because they assumed isothermal atmosphere 
without a steep density gradient at the disk-halo interface, and partly because of the lower 
numerical resolution. Since the downflow speed exceeds the local sound speed, strong shock 
waves are formed at the magnetic loop footpoints. Figure 5b shows that the shock waves heat 
the cool gas around z/r ~ 0.1. 

At the final phase (t/to — 1-6), the expansion of the magnetic loops is stalled because 
the driving forces diminish at the top of the loop. Figure 6 shows the vertical distribution 
of the magnetic pressure, gas pressure, (3, and density at the midpoint of the emerging loop 
(x/ro = 0.25). From this figure, it is clear that the magnetic pressure and gas pressure gradients 
nearly balance at the top of the emerging loop. This magneto-hydrostatic state is attained 
because the Parker instability is stabilized in the hot layer where the local pressure scale height 
becomes larger than the length of the magnetic loop. 

Figure 7 shows the snapshots of the density distribution for all models (/3q = 0.5,1,2,4 
and 10) at the stage when the top of the magnetic loops enters the hot region (z/r > 0.24). 
The time scale for the loop emergence is shorter for lower (3, consistent with the results of the 
linear analysis presented in section 2.3. Figure 8 shows the vertical distributions of the density 
at the midpoint of the emerging loop. The density increases with height at the top of the 
magnetic loops in all models. In the model with (3q = 10, the emerging loop stalls at a lower 
height (z/r < 0.3) because the released magnetic energy is small. 

4. Discussion 

4-1. Formation of Dense Loop Structure 

Let us discuss why the emerging loop becomes top-heavy. Figure 9 shows the magnetic 
field lines at t/to = 1.4 (solid curves) and t/to = 0.9 (broken curves) for the fiducial model (model 
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B). These field lines are iso-contours of the y-component of the vector potential. The magnetic 
field lines comoving with the plasma in ideal MHD can be identified by the value of the vector 
potential. Since the frozen-in condition is assumed and the numerical diffusion is negligibly 
small, the field lines at t/to = 0.9 (broken curves) have moved with gas to the corresponding 
field lines at t/to — 1-4 (solid curves). As the magnetic loops rise, they can accumulate the gas 
above the loop when the loop top is flat. Figure 9 shows that the loop shape at t/to = 0.9 is 
favorable for the mass accumulation around the loop top. At t/to = 1-4, since the curvature 
increases, the mass accumulated around the loop top slides down along the magnetic field lines. 
Flat top loops are formed when the loop length is much longer than the local pressure scale 
height. In our model atmosphere, since the local pressure scale height decreases with height 
in < z/r < 0.1 and increases with height in 0.1 < z/r < 0.25 (see Figure 2b), the loop top 
tends to be flat when its height is z/ro ~ 0.1. The dense shell around the loop top is formed 
when 0.1 < z/ro < 0.2, and the mass drains as the loop rises. 

In order to quantitatively evaluate the density enhancement at the top of the emerging 
loop, we adopt the same method of analysis as Isobe et al. (2006). The velocity vector v can 
be divided into two components: 

V = V± + V\\, (11) 

where v± and 

*=w\ B ' < 12 > 

are the velocity components perpendicular and parallel to the magnetic field, respectively. 

Figure 10 shows the distribution of v and V ■ v (Figure 10a,b), v± and V • (Figure 
10c, d), and v» and V • v» (Figure 10e,f) at t/to — 1-2 and t/to — 1-3 for the fiducial model. The 
distributions of V ■ Vu and V ■ v± show that the density inside the emerging magnetic loops 
keeps decreasing. However, V ■ v and V • v± at the top of the emerging loop are negative. 
Therefore, at least for the present parameters, it can be concluded that the density increases 
in the top region of the emerging loop. Figure 11 shows the distribution of the density near 
the loop top at t/to = 11 and t/to — 1-2 for the fiducial model. The clump of gas near the loop 
top survived due to the small curvature in the low-temperature layer and is later compressed 
in the high-temperature layer. 

4-2. Comparison with Observations 

Fukui et al. (2006) found two dense gas features having a loop-like shape with a length 
of several hundred pc within ~ 1 kpc from the Galactic center. Figure 3a shows that the 
most unstable wavelength is likely to be 350 - 500 pc. Our numerical simulation produced the 
magnetic loops whose wavelength is ~ 400 pc and their height is ~ 350 pc. In the observed 
molecular loops, the line-of-sight velocity along the loop changes linearly with the arc-length 
and has large gradients (~ 30 km s -1 ). Figure 12 shows the velocity components along the 
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outermost magnetic field line at t/to = 1.4 for the fiducial model. The gradient of the downflow 
velocity given in Figure 12c is ~ 20 — 30 km s _1 . Since the downflow speed exceeds the local 
sound speed, strong shock waves are formed at the magnetic loop footpoints. Such shocks can 
be the origin of the observed large velocity dispersions near the footpoints of Galactic center 
molecular loops. 

4-3. Estimation of the Kinetic Energy of the Emerging Loop 

Fukui et al. (2006) estimated that the kinetic energy of the molecular loop is ~ 10 51 erg. 
In our simulation, the kinetic energy flux, F^, carried by the downflow is 

F k ~ -Pdfvlt ~ 10~ 5 erg cm' 2 b _1 , (13) 

where pdf is assumed to be 10~ 24 g cm -3 and v^f is found to be ~ 30 km s _1 according to the 
numerical result for the fiducial model. The numerical results indicate that the loop width in 
z-direction is about 30 pc (~ 10 20 cm) and the downflow continues for about 10 15 s ~ 3 x 10 7 
yr. When the loop thickness in y-direction is 100 pc, the total kinetic energy of the downflow 
toward the both footpoints is ~ 10 51 erg. This energy is comparable to that estimated from 
the observation. 

4-4- Effects of Cooling and Rotation 

In this study, we assumed that the gas layer at the initial state is composed of the cool 
equatorial layer, warm layer, and hot layers. We assumed adiabatic gas but in the interstellar 
gas, cooling and heating play essential roles in the formation of molecular clouds (Field 1965). 
Kosihski & Hanasz (2006, 2007) investigated the Parker instability coupled with thermal pro- 
cesses (cooling and heating). They found that the Parker instability can trigger the thermal 
instability which form dense clouds in the valleys of the magnetic loops. In subsequent papers, 
we would like to include the gas cooling and heating effects. 

We should consider the effect of the rotation of the disk to construct a more realistic 
model. In rotating disks, the Parker instability is slightly suppressed by Coriolis forces (Chou et 
al. 1997). Hanasz et al. (2002) presented the results of resistive 3D MHD simulations of a local 
part of the disk including the contribution of the Coriolis force. They demonstrated that the 
Parker instability, the twisting of the loops, and magnetic reconnection lead to the formation 
of helically twisted magnetic flux tubes. It is worth noting other 3D effects such as the growth 
of the interchange mode (e.g., Nozawa 2005). Isobe et al. (2006) showed that dense shells at 
the top of magnetic loops fragment into filaments due to the growth of the Rayleigh- Taylor 
instability. In order to include this effect, we need to carry out 3D MHD simulations. 

In Galactic gas disks, magneto-rotational instability (MRI; Balbus & Hawley 1991) grows 
and drives magnetic turbulence inside the disk. In the nonlinear stage, MRI drives magnetic 
turbulence inside the disk. The effects of stochastic magnetic fields on the growth of the Parker 
instability was reported by Parker & Jokipii (2000). The results of global 3D MHD simulations 
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of the galactic center gas disk are reported by Machida et al. (2009). 

4-5. Dependence on the initial disk temperature 

We studied the dependence of the numerical results on the initial disk temperature. 
Figure 13 shows the density distribution at the stage when the top of the magnetic loops 
reaches Z/r ~ 0.35 for models with lower initial disk temperature (T c = 0.05T ) and higher 
initial disk temperature (T c = 0.2T ). Other paremeters are the same as the fiducial model 
(model B). Numerical results indicate that the equatorial dense region at this stage is thin for 
the cool disk (T c = 0.05T ) and thick for the warm disk (T c = 0.2T ). Although the length of 
the loops near the equator is smaller for cooler disk, the length of the loops emerging in the 
halo is almost the same (~ 400 pc). This is because the most unstable wavelength of the Parker 
instability is determined by the local pressure scale height. 

4.6. Other Effects 

Finally, we briefly mention about other effects such as cosmic rays, self-gravity and the 
presence of a spiral arm of the Galactic gas disk. In Galactic disks, since the cosmic ray pressure 
is comparable to the gas pressure, cosmic rays enhance the growth rate of the Parker instability 
(e.g., Parker 1966, Hanasz & Lesch 2000, 2003; Kuwabara, Nakamura & Ko 2004). The effects 
of the cosmic rays may be more important in the Galactic center where strong activities can 
produce high energy particles. 

In this paper, we neglected the self-gravity of the gas. When the surface density of 
the gas disk is large enough, Parker- Jeans instability grows (e.g., Elmegreen 1982; Nakamura, 
Hanawa & Nakano 1991; Kim, Ostriker & Stone 2002; Lee et al. 2004). The Jeans instability 
has a larger growth rate than does the Parker instability when either the magnetic field is weak 
or the wavelength of the perturbation is long. 

Franco et al. (2002) showed by 3D MHD simulations that the Parker instability creates 
massive clouds inside the spiral arm of the Galactic gas disk and that the dense gas accumulated 
around the equatorial plane forms corrugated structure. The distribution of H I gas and dust 
below the molecular loops found in the Galactic center (Torii et al. 2009) is consistent with 
the simulation by Franco et al. (2002). The density distribution in the nonlinear stage of our 
simulation also shows such corrugated structures. 

5. Summary 

In this paper, we have carried out 2D MHD simulations of the Galactic center gas disk 
consisting of the low-temperature layer and the mid-temperature layer with the overlying hot 
halo. We found that the numerical results reproduce basic features in the Galactic-center 
molecular loops as observed with NANTEN, such as the loop length, the velocity gradient 
along the loops, and large velocity dispersions at the footpoints of the loops. 
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We also discussed the reason why the top of the emerging loop becomes over-dense. This 
is because the effective gravity along the magnetic field lines decreases when the curvature radius 
of the magnetic field lines increases in the low-temperature layer. The gas survived near the 
loop top due to its small curvature is compressed when the loop top enters the high-temperature 
layer. 

Finally, we suggest that the Galactic center molecular loops are analogous to the arch 
filaments in the solar chromosphere. The molecular loops emerge from the low temperature 
layer just like the dark filaments observed in Ha images of the emerging flux region of the sun. 

We are grateful to T. Sakurai, T. Kudoh and D. Shiota for useful comments and dis- 
cussion. Numerical computations were carried out on the general-purpose PC farm at Center 
for Computational Astrophysics, CfCA, of National Astronomical Observatory of Japan (P.I. 
KT). This work is financially supported in part by a Grant-in- Aid for Scientific Research 
(KAKENHI) from JSPS (No. 20244014). This work is also carried out by the joint research 
program of the Solar- Terrestrial Environment Laboratory, Nagoya University. 
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Table 1. Units for Normalization 



Physical Quantity Symbol Normalization Unit 



Physical Value 



Length x,z r 

Density p p 

Velocity v Cso 

Pressure P PqCsi 



3 x 10" 10 



1 kpc 
1CT 22 g cm~ 3 



18 km s 



dyne cm 



Magnetic Field B (A)Cf ) 1/2 

Temperature T fiCg /(^R) 



17 n G 
2 x 10 4 K 
5.6 xlO 7 yr 
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Fig. 1. Schematic picture of the simulation model and simulation box. 
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Fig. 2. (a) Distribution of the initial density (dashed curve), gas pressure (dotted curve), magnetic pres- 
sure (solid curve), and temperature (dash dotted curve), for the fiducial model (model B). (b) Distribution 
of the initial local pressure scale height A (solid curve) and modified scale height Ab = (l+/3(z) _1 )A (dashed 
curve) for the fiducial model (model B). 
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Table 2. List of models 



Model 


r [kpc] 


zi/r 


22 / r 


Zf/ro 


Po 


A 


1 


0.12 


0.24 


0.08 


0.5 


B a 


1 


0.12 


0.24 


0.08 


1 


C 


1 


0.12 


0.24 


0.08 


2 


D 


1 


0.12 


0.24 


0.08 


4 


E 


1 


0.12 


0.24 


0.08 


10 



Fiducial model 
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Fig. 3. (a) Linear growth rate of the Parker instability for the unperturbed states with (3q = 0.5,1,2,4 
and 10 as a function of wave number fe x . (b) The dependence of the maximum growth rates on the initial 

— 1/2 

plasma beta. The solid line shows a line where iu> oc /? . The unit of the growth rate is CsoA'o = l/^o- 
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Fig. 4. Numerical results for the fiducial model (model B); density in a logarthmic scale (colors), magnetic 
field lines (solid curves), and velocity field (vectors) for t/t = 0.0,0.5, f.0, f. 2, 1 .4 and 1.6. 
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Fig. 5. Numerical result for the fiducial model (model B; /3q = 1); (a) density in a logarithmic scale 
(colors) at t/to= 1.4. (b) temperature in a logarthmic scale (colors), magnetic field lines (solid curves), 
and velocity field (vectors) at t/t = 1.4. 



16 



t/t =1.6 |S =1 X/r =0.25 
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Fig. 6. Onc-dimensional (z-) distribution of the logarithmic magnetic pressure P m (solid curve), gas pres- 
sure P g (broken curve), plasma beta (3 (dotted curve) and density p (dash dotted curve) at t/to = 1.6 for 
the fiducial model (model B). 
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Fig. 7. Numerical results for all models {(3q = 0.5,1,2,4 and 10); density in a logarthmic scale (colors), 
magnetic field lines (solid curves), and velocity field vectors. 
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Fig. 8. Distribution of the logarithmic density in the z-direction measured at the loop center for models 
(a) model A, (b) model B, (c) model C, (d) model D and (e) model E at various times. 
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Fig. 10. (a,b) Divergence of velocity at t/to = 1.2 and 1.3. Solid curves show magnetic field lines, and 
arrows show the velocity vector, (c,d) for v±, and (e,f) for v\\. Colors in figures (c,d) and (e,f) show the 
divergence of velocity components perpendicular and parallel to the magnetic field, respectively. 
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Fig. 11. Distribution of the logarithmic density (solid curve) and pressure (broken curve) near the loop 
top at x/tq = 0.25 for the fiducial model (model B) at t/to = 1.1 (left) and t/to = 1.2 (right). 
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Fig. 12. (a) Numerical results for the fiducial model (model B); density in a logarthmic scale (colors), 
magnetic field lines (solid curves), and velocity field vectors at t/to = 1 .4. (b) The outermost magnetic field 
line at t/to = 1-4. (c) Velocity components (v x ,v z ) along the outermost magnetic field line, (d) Density 
along the outermost magnetic field line at t/to = 1.4 (solid curve) and t/to = (broken curve). 
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Fig. 13. Dependence of the numerical results on the initial disk temperature. Color shows the density 
distribution at the stage when the top of the magnetic loops reaches z/ro ~ 0.35. (a) T c = 0.05Tb at 
t/t = 1.0. (b) T c = 0.2T a t t/t = 2.0. Solid curves show magnetic field lines. Arrows show velocity 
vectors. 
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